Gas Doping on the Topological Insulator Bi2Se3 Surface 
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Gas molecule doping on the topological insulator Bi2Se3 surface with existing Se vacancies is 
investigated using first-principles calculations. Consistent with experiments, NO2 and O2 are found 
to occupy the Se vacancy sites, remove vacancy-doped electrons and restore the band structure of a 
perfect surface. In contrast, NO and H2 do not favor passivation of such vacancies. Interestingly we 
have revealed a NO2 dissociation process that can well explain the speculative introduced "photon- 
doping" effect reported by recent experiments. Experimental strategies to validate this mechanism 
are presented. The choice and the effect of different passivators are discussed. This step paves the 
way for the usage of such materials in device applications utilizing robust topological surface states. 

PACS numbers: 71.20.-b,73.20.-r,71.20.Nr 



Three-dimensional (3D) topological insulators (TI) 
have attracted extensive research interest recently \^\^. 
Their novel topological surface states (TSS) in the bulk 
energy gap open great practical potential in spintron- 
ics and quantum computation by realizing the Majorana 
fermions [5J. However, most TI materials available today 
are poor insulators in the bulk due to heavy defect dop- 
ing jsl, 13, and references therein] , hindering the utiliza- 
tion of the TSS in practice. To remove the bulk carriers, 
surface doping with gas molecules was adopted as a pow- 
erful tool in experiments though the underlying 
detailed atomistic mechanisms are still unknown. Ex- 
ploring chemical reactions and understanding the modi- 
fied surface structures are crucial to comprehend the cur- 
rent experimental results and help finding new methods 
to achieve the bulk insulating state. 

As of todw, the most attractive TI material is 
Bi2Se3 @, [nl, demonstrating simple Dirac-type TSS 
along with large bulk energy gap. Presumably due to the 
Se vacancies, however, it is found to be a n-doped semi- 
conductor. To achieve a real bulk insulator, NO 2 [l,[E[l3l 
and O2 0, [3 g^'S species were used to dope the surface 
and successively remove the donated bulk conduction 
electrons. In particular, the NO2 doped Bi2Se3 surface 
exhibits a mysterious behavior, in which the surface loses 
electrons when exposed to photon flux in angle-resolved 
photoemission spectroscopy (ARPES) experiments [1,J3- 
This is very different from the cases of graphene [l3| 



and another known TI material Bi2Te3 [8|, |14| |. where 
the surface gains electrons back through desorbing NO2 
by the photon exposure. Though this photon-assisted 
stimulation or the so-called "photon-doping" method was 
employed to manipulate the surface bands in experi- 
ments 15[, the underlying mechanism remains to be de- 



fined. In this letter, we have studied from first-principles 
the reaction of gas molecules on the Bi2Se3 surface con- 
taining Se vacancies and conclude how these affect the 
TSS. 

The first-principles molecular dynamic (FPMD) simu- 




FIG. 1. (a) Schematic of NO2 molecule adsorption and dis- 
sociation processes and (b) corresponding reaction energetics 
on the Bi2Se3 surface with a Se vacancy. The Se vacancy 
is indicated by the white dot (step I) and donates two elec- 
trons that can dope the surface. During NO2 exposure, the 
molecule occupies the vacancy site quickly with one O binding 
to three Bi atoms from the second atomic layer (step II). The 
corresponding N-O bond become weaker compared to before 
adsorption (step I). The stable adsorption structure is shown 
in (b) together with charge density difference before and af- 
ter adsorption. The donor electrons transfer partially from 
the vacancy to NO2. Under external photon exposure (step 
III), the weakened N-O bond breaks, resulting in the NO2 
dissociation into NO -I- O (step IV). The NO molecule leaves 
the surface, while O passivates the vacancy. In this case, two 
donor electrons transfer to the O atom and the vacancy is 
compensated. White balls stand for Bi atoms, red for O, and 
blue for N. In the isovalue surface plots of the charge differ- 
ence, solid blue (wired gray) color denotes charge depletion 
(gain). 



lation was employed to investigate the dynamic process 
of chemical adsorption and search for the stable atomic 
configurations. In bulk Bi2Se3, five atomic layers form 
a quintuple layer (QL) while the coupling between two 
such QLs is of the van der Waals type [11]. Here we 
used a two-QL-thick slab model containing 3x3 primi- 
tive unit cells in the xy plane to simulate a surface for 
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FIG. 2. Atomic structures and band structures for (a) a perfect Bi2Se3 surface, (b) a surface with a Se vacancy, (c) NO2 
adsorption, (d) O adsorption from the dissociated NO2, (e) O2 and (f) H2 adsorption. Bands are aligned to their valence band 
maximum at F. The Fermi energy is shifted to zero.(a)-(d) show structures from side view, whereas (e)-(f) from the top view. 
Red lines indicate the surface states. White balls stand for Bi, yellow for Se, red for O, blue for N and small gray for H atoms. 
The transparent ball indicate the missing Se at the vacancy in (b) . Band structures are calculated from a 3 x 3 supercell of 
seven QLs thick. Colored circles denote Dirac points. 



the FPMD adsorption studies. Thicker slabs were tested 
to yield correct and converged results. Both the pristine 
surface and the surface with a Se vacancy at the out- 
most layer were considered for adsorption of NO2, NO, 
O2 and H2 molecules. FPMD calculations including van 
der Waals interactions [l^ have been performed at 300 K 
in microcanonical ensemble (NVE) with an integration 



time-step of 1.0 fs [17|. We let the calculations run for 



5.0 ps, while in worst case, roughly after 1.5 ps no more 
significant structural changes were observed. Using the 
optimized structures from FPMD, the slab was extended 
to seven QLs thick to prevent the interactions between 
states from the top and the bottom surfaces [l3, [l9| . 
Inversion symmetry was maintained by placing the same 
molecules on both surfaces. Subsequently, band structure 
calculations were carried out to investigate the electronic 
properties. (2^ 

The results are shown in Fig. [21 The pristine surface 
has a single pair of Dirac-like TSS with the Fermi level 
Ep crossing the Dirac point at F, consistent with pre- 
vious calculations 
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When a Se vacancy forms on 
the surface, Ep shifts to the bulk conduction band bot- 
tom, resulting in n-type doping as observed in experi- 
ments. Interestingly, the surface states are energetically 
changed considerably by the existing vacancies. The orig- 
inal Dirac point at F shifts upward and almost merges 
into the conduction band, while a new Dirac point forms 
at the M point (Fig. E^b)). In the bulk gap, Ep still 
intersects the TSS odd times, indicating the topologi- 
cal nontrivial character. It should be noted that in our 
model Se vacancies distribute periodically on the surface. 



reflecting the employed supercell periodic boundary con- 
dition. However due to the vacancy random distribu- 
tion on the real sample surfaces, the ARPES measured 
only a single Dirac point at F in the primitive Brillouin 



zone [y, llJ, |lj|. This illustrates a possible way to en- 
gineer the Dirac point by designing well-ordered surface 
potential. 

NO2 and NO. In the FPMD simulations, NO2 
molecules are placed above both the pristine surface and 
the surface with vacancy. NO2 moves freely above the 
pristine surface whereby no strong binding effect is ob- 
served. The total energy calculation reveals that NO2 
physisorbs on this surface with a weak adsorption poten- 
tial energy (APE) of - 0.4 eV, resulting in tiny charge 
transfer from the surface to the molecule. However, on 
the defected surface NO2 chemisorbed at the vacancy site 
quickly. One of the NO2 oxygens, labeled as 01, adsorbs 
and binds to three Bi atoms from the second atomic layer, 
see Fig.[lja) step II. The equilibrium bond length of N-01 
is elongated to 1.6 A compared to the other N-0 bond of 
1.3 A, indicating weakened N-01 bond upon adsorption. 
Proceeding further, different FPMD replicas resulted into 
two situations. Firstly, in the majority of the simulations, 
stable adsorption of NO2 at the vacancy with the APE of 
- 2.0 eV took place. The corresponding band structure is 
depicted in Fig.[2Kc), whereby the Dirac cone is recovered 
at F. However, the surface is yet slightly electron-doped 
with Ep 0.15 eV above the Dirac point. This indicates 
that NO2 can not take away all donor electrons from 
the vacancy. In the other situation, localized vibrations 
cause the weakened N-01 bond to break. Subsequently, 
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NO2 dissociates into NO which leaves the surface, while 
the 01 atom passivates the vacancy. The respective APE 
decreases to -1.6 eV. To estimate the dissociation energy 
barrier, we performed density-functional calculations us- 
ing the nudged elastic band method 2l[ and monitor 



dissociation to happen smoothly at a barrier of about 
0.4 eV. From the related band structure in Fig. [2jd), we 
see the 01 atom passivation shifts Ep further down to 
the Dirac point and totally restores the band structure of 
the pristine surface. This can be understood by the fact 
that the single 01 atom can accommodate all two extra 
electrons from the Se vacancy, more than NO 2 in the for- 
mer case, as illustrated by the charge density change in 
Fig. mb). Contrary to NO2, NO binds weakly (APE of 
-0.6 eV) and does not yield a similar good passivation of 
the vacancy, due to its much lower electron affinity. 

It is interesting to compare these results with recent 
ARPES measurements In experiments, the Fermi 

level is 0.3 eV above the Dirac point. Supposing all doped 
electrons are due to the Se vacancies, a dosage of 0.1 
Langmuir(L) of NO2 can shift Ep down by 0.15 eV. Af- 
ter dosing with NO2, in experiment two additional non- 
dispersive peaks appear at the binding energies of -4.0 eV 
and -7.5 eV Q. This agrees with our calculated density 
of states (DOS) projected onto NO2 (Fig. [HJb)), where 
partial DOS analysis (not shown here) indicates that the 
peak at ^ -4.0 eV is mainly due to both O atoms and the 
other at ~ -7.5 eV due to all N and O atoms. The ARPES 
results are found to remain stable with temperature vary- 
ing from 10 K to 300 K. This can be explained by the 
large APE and the dissociation barrier (0.4 eV), much 
too high to be overcome by RT thermal fluctuations. 
Thereby NO2 will neither desorb nor dissociate at 300 K. 
Moreover, the surface with a dose of 0.1 L can also be fur- 
ther hole-doped by extensive photon (energy 28-55 eV) 
exposure until Ep reaches the Dirac point [15|, called as 
"photon-doping" . Considering that NO2 molecules disso- 
ciate into NO -I- O under laser pulse irradiation |22] , the 
above experimental observation may be explained by the 
NO2 dissociation process (Fig. [1] step III). High-energy 
photons stimulate the N-01 bond to break and subse- 
quently transfer remaining electrons from the vacancy to 
the 01 atom, finally downshifting Ep to the Dirac point. 
We conclude that a dosage of 0.1 L of NO2 passivates 
most of the Se vacancies on the surface, as the experi- 
mental Ep position (0.15 eV above the Dirac point) is 
the same as our calculated one, see Fig. EJc). Therefore 
photon exposure induces an 0-passivated surface without 
extra donor electrons. In addition, without photon expo- 
sure a dosage of 2.0 L more can also shifts Ep down to 
the Dirac point. This can be attributed to the weak ad- 
sorption of NO2 on the pristine surface, see above, where 
much less electron transfer happens compared to that of 
the strong adsorption on the vacancy. 

The proposed dissociation mechanism can be easily 
checked in experiments. One way is to monitor whether 




FIG. 3. Density of states (DOS) for (a) the perfect pristine 
surface and the one with a Se vacancy, (b) NO2, (c) a single 
O atom, (d) O2, and (e) H2 passivated. To be more clear the 
adsorbates PDOS have been scales by Nx/Nto*, where Nx is 
the number of atoms in the adsorbate and Ntoi is the total 
number of atoms in the supercell. 



the composition of NO (NO2) increases (decreases) af- 
ter photon exposure. The second way is to measure the 
valence band spectra. We predicted that the peak at 
~ -7.5 eV will become weaker or even disappear after ex- 
tensive photon exposure, as indicated by the projected 
DOS in Fig. (He). The third alternative is to see the sur- 
face using STM or from vibration spectra. Additionally, 
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one can find that the monoatomic oxygen, which is gen- 
erally used in semiconductor industry for plasma ashing, 
can be an excellent passivator to Se vacancies, if under 
controllable exposure. 

O2 and H2- The O2 molecule also gets adsorbed in the 
vacancy site quickly with an estimated APE of - 3.4 eV. 
The two O atoms still remain connectted to each other 
with one O (labeled as 01) binding to one Bi atom and 
the other O (02) binding to the other two Bi atoms, as 
Fig. [2Ke) shows. In the equilibrium structure, the 01-Bi 
bond length is ^ 2.2 A, while two 02-Bi bonds form at 
~ 2.4 A. Concluding from the band structure analysis, 
O2 captures two extra electrons from the vacancy and 
recover the neutral-charged TSS very well. After adsorp- 
tion no O2 dissociation was observed in any of our MD 
simulations, clearly due to stiffness of 0=0 bond. This 
suggests why photon exposure had not any improving ef- 
fect in experiments with O2 0, @1- In contrast, H2 does 
not occupy the vacancy site in the MD adsorption sim- 
ulations. However after artificially saturating dangling 
bonds with hydrogens, the relaxed structure shows simi- 
lar bondings with three below Bi atoms as in the O2 case, 
see Fig. [Iff). This passivation can also restore a perfect 
Dirac-type band structure (Fig. [Sff)). Though H2 ad- 
sorption is not favored in dynamics, single H atoms can 
be an ideal passivation to Bi dangling bonds in theoret- 
ical studies, removing nontrivial surface dangling bond 
states and simplify the physics underlying. 

In summary, we have studied the chemical adsorption 
of gas molecules (NO2, NO, O2 and II2) on the Bi2Se3 
surface. The passivation of common Se surface vacan- 
cies is found to gain the donor electrons much more ef- 
fectively than weak adsorption on the pristine surface. 
The local atomic structures after adsorption are revealed 
and the band structures are compared. NO2 is observed 
to passivate the vacancy and accommodate partially the 
vacancy-doped electrons. By overcoming a moderate en- 
ergy barrier {e.g. with laser stimulation), NO2 dissoci- 
ates into NO and a single O atom, whereby the latter 
occupies the vacancy and accommodates all donor elec- 
trons. This dissociation and charge transfer can explain 
the mysterious "photon-doping" effect seen in ARPES 
experiments. As well, both O2 molecule and the single 
O atom are found to passivate the vacancy very well and 
restore the surface band structure to the charge neutral 
state. In contrast, NO and II2 are not favorable adsor- 
bates on the Bi2Se3 surface. The combinination of MD 
with electronic structure analysis is shown to be versatile 
in explaining experimental observations along with sug- 
gesting new routes for device engineering, by tailoring 
different dopants and TI materials. 
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